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Abstract: Determinations of structure functions and parton distribution functions have 
been recently obtained using Monte Carlo methods and neural networks as universal, un- 
biased interpolants for the unknown functional dependence. In this work the same meth- 
ods are applied to obtain a parametrization of polarized Deep Inelastic Scattering (DIS) 
structure functions. The Monte Carlo approach provides a bias-free determination of the 
probability measure in the space of structure functions, while retaining all the information 
on experimental errors and correlations. In particular the error on the data is propa- 
gated into an error on the structure functions that has a clear statistical meaning. We 
present the application of this method to the parametrization from polarized DIS data 
of the photon asymmetries A\ and Af from which we determine the structure functions 
gP(x,Q 2 ) and gf(x,Q 2 ), and discuss the possibility to extract physical parameters from 
these parametrizations. This work can be used as a starting point for the determination 
of polarized parton distributions. 
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1. Introduction 

In QCD the description of scattering processes at large momentum transfer (Q 2 ^> Aq^) 
involving (polarized) hadrons in the initial state is based on the factorization theorem. 
The latter allows a separation between the high-energy dynamics, described by coefficient 
functions which are calculable in perturbative QCD, from low-energy, non-perturbative ef- 
fects, binding partons into hadrons, which are encoded into (polarized) parton distribution 
functions (PDFs). 

The growth in statistics and increase in precision of data from experiments involving 
polarized hadrons scattering calls for a more accurate determination of polarized PDFs 
and their errors. A crucial problem in this respect is the determination of the uncer- 
tainty on a function (i.e. a probability measure on a space of functions) from a finite set 
of experimental data points. In the standard PDF extraction approach to the problem 
the infinite-dimensional space of continuous functions is mapped into a finite-dimensional 
space of parameters by choosing a particular basis in the space of functions and truncating 
the basis to a finite number of elements. This procedure entails some degree of arbitrari- 
ness. Any sensible choice must strike a balance between two competing requirements: on 
the one hand a small number of parameters introduces a bias in the determination of both 
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the functional form and the errors, as the chosen parametrization would not allow enough 
flexibility; on the other hand a large number of parameters could spoil the convergence of 
the fit, or be too sensitive to the statistical fluctuation of the experimental data. 

This problem has been addressed by the NNPDF Collaboration in the case of unpo- 
larized Deep Inelastic Scattering (DIS) structure functions in Refs. [1,2], and in the case 
of the unpolarized PDFs in Refs. [3-5] using a method based on statistical inference and 
neural networks as an interpolating tool. 

While avoiding technical complications linked to the extraction of PDFs from ob- 
servables, the determination of structure functions addresses the main issue of devising a 
faithful estimation of errors on a function extracted from experimental data. The main 
ingredient in the studies above is the usage of Monte Carlo methods to obtain a repre- 
sentation of the probability measure in the space of structure functions. An ensemble of 
artificial data is generated, which reproduces all the statistical features (i.e. variances and 
correlations) of the original experimental data. Each set of artificial data is called a replica. 
A structure function, parametrized by a neural network, is then fitted to each replica. The 
net result of this procedure is an ensemble of fitted functions. This ensemble of fitted func- 
tions provides a representation of the measure in the space of structure functions. Errors 
and correlations of any observable involving the structure functions are obtained averaging 
over the ensemble of fits. Moreover suitable statistical estimators can be defined from the 
Monte Carlo ensemble which provide a quantitative description of the possible biases and 
inconsistencies in the fitting procedure. This method has been described in great detail in 
Refs. [1-4] to which the interested reader should refer. 

The aim of this work is to apply the same techniques to obtain a bias-free parametriza- 
tion of the photon asymmetries A\ and Af from available polarized DIS data and extract 
from them the corresponding structure functions g\ and gf . We provide further testing of 
the Monte Carlo method, and produce statistically meaningful error bars for the structure 
function. Besides allowing us to address all systematics related to the data and the method, 
such a parametrization might be an ideal input for a fit based on factorization scheme- 
invariant evolution equations to determine a s , as proposed in Refs. [6,7]. As shown in this 
work, a careful treatment of statistical and systematic errors leads to a reliable extraction 
of physically meaningful parameters such as a s , qa , and the higher-twist contributions to 
the structure functions. While these are not the best determinations available for these 
parameters, the results we obtain are in agreement with other determinations, and show 
the robustness of the Monte Carlo method. 

We shall now discuss in turn the two steps that are needed to produce the Monte 
Carlo sample of fitted functions: first the treatment of the experimental data, and then the 
actual fitting procedure. The experimental data points included in the fit are discussed in 
Section 2; Section 3 summarizes briefly the NNPDF approach and the characteristics of 
the neural networks used for this particular study. The results of our fits, together with 
their phenomenological implications are presented and discussed in Sections 4 and 5. 
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2. Experimental data 



The cross section asymmetry for parallel and anti-parallel configurations of longitudinal 
beam and target polarizations is given by: 

An = -r. =r 2.1 

and it is related to the virtual-photon asymmetries A\ , A2 by: 

Ay =D(Ax+ V A 2 ) ~DA X . (2.2) 
The photon depolarization factor D depends on kinematic factors and on the ratio: 



where ai and <jt are the longitudinal and the transverse cross sections respectively (see 
e.g. Refs. [8,9] for a detailed definition of all the quantities). 

The polarized structure functions g\ and g 2 are related to the virtual-photon asymme- 
tries by: 

A 2 {xM)=l , (2.5) 

where 

F ^ = Mtr£m F2{x - Q2) - (2 ' 6) 

is the unpolarized structure function, 7 2 = in ^-f , and m denotes the nucleon mass. 

The main features of each experimental data set used in the present analysis are 
summarized in Tab. [fl, and their kinematical coverage of the (x, Q 2 )-plane is shown in 
Fig. |]. We observe that the kinematical coverage of the available data is rather small, 
especially when compared to the one of the available unpolarized DIS data, thus we will 
have a sizable region of the kinematical plane in which the fit extrapolates the behaviour 
extracted from the region covered by data. 

From Tab. || we infer that the systematic errors are on average one order of magnitude 
smaller than the statistical ones. This justifies the procedure of neglecting correlations for 
systematic errors and the procedure of summing errors in quadrature when computing the 
figure of merit (x 2 ) to be minimized in the fitting procedure. 

Finally we notice that El 55 data have been corrected to yield A± by adding in Eq. 



(2.4) the 52 contribution evaluated with the Wandzura-Wilczek relation and using the 
parametrization of g\/F\ given in Ref. [14]: this shift is also added as a source of uncertainty 
in the total error of the data set. 
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Experiment 


x range 


Q 2 (GeV 2 ) range 




(ffstat) 


(°"syst) 


{o"norm ) 


Type 


Ref. 


Proton 


EMC 


0.015 - 0.466 


3.5 - 29.5 


10 


0.077 


0.024 


0.028 


Ax 


[10] 


SMC 


0.001 - 0.480 


0.3 - 58.0 


15 


0.026 


0.003 


0.012 


Ax 


[11] 


SMC low-z 


0.0001 - 0.121 


0.02 - 23.1 


15 


0.033 


0.002 


0.006 


Ax 


[12] 


E143 


0.031 - 0.75 


1.27 - 9.52 


28 


0.045 


0.016 


0.012 


Ax 


[13] 


E155 


0.015 - 0.75 


1.22 - 34.72 


24 


0.043 


0.018 


0.026 


9i /Fx 


[14] 


HERMES06 


0.0058 - 0.7311 


0.26 - 14.29 


45 


0.126 


0.019 


0.017 


Ax 


[15] 


Deuteron 


COMPASS 


0.0051 - 0.474 


1.18 - 47.5 


12 


0.034 


0.017 


0.011 


Ax 


[16] 


SMC 


0.001 - 0.480 


0.3 - 58.0 


15 


0.032 


0.003 


0.006 


Ax 


[11] 


SMC low-z 


0.0001 - 0.121 


0.02 - 23.1 


15 


0.069 


0.005 


0.005 


Ax 


[12] 


E143 


0.031 - 0.75 


1.27 - 9.52 


28 


0.066 


0.011 


0.008 


Ax 


[13] 


E155 


0.015 - 0.75 


1.22 - 34.72 


24 


0.091 


0.009 


0.011 


9i /Fx 


[14] 


HERMES06 


0.0058 - 0.7311 


0.26 - 14.29 


45 


0.089 


0.007 


0.009 


Ax 


[15] 



Table 1: The proton and deuteron experimental data sets included in the present analysis. We show 
the kinematic range, the number of points, the average statistical, systematic and normalization 
uncertainty, and the measured observable. 
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Figure 1: Experimental data in the (x, Q 2 ) plane used in the present analysis for the proton (left) 
and for the deuteron (right) target. 



3. The NNPDF approach 



In this section we briefly review the approach used to extract an unbiased determination of 
the asymmetry A± and the structure function g\ from the available inclusive polarized DIS 
data, following the analysis performed by the NNPDF Collaboration for the determination 
of the unpolarized structure function [1,2] and the parton densities [3-5]. 

The core idea underlying the NNPDF approach is based on using Monte Carlo methods 
to build a representation of the probability measure in the space of structure functions, 
and parametrizing the space of structure functions using neural networks. We refer the 
interested reader to the papers cited above for a detailed description of the methods and 
in the following we will briefly discuss the settings used in this analysis. It is worthwhile 
emphasizing that the Monte Carlo method does not require the use of neural networks, 
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and would yield a robust determination of the errors with any parametrization of the 
structure function, provided the parametrization is sufficiently flexible. A comparison of 
Monte Carlo analyses based on different parametrizations was performed in the framework 
of the HERA-LHC workshop by comparing the standard HI and NNPDF analyses. The 
greater flexibility of the neural network parametrization compared to fixed functional forms 
is reflected in the larger error bands obtained using the NNPDF method, especially when 
considering the x region not covered by data (i.e. the extrapolation region). These features 
are illustrated by the results in Sects. 3.2 and 3.4 of Ref. [17]. 

3.1 Monte-Carlo replicas 

We generate N Tep Monte-Carlo replicas of the experimental data according to 

A<r ) ' k (x,C?)=(l + r kfN ^ ) [At P) (x,Q 2 ) + r k ,M*,Q% (3-1) 
V A 1 (x,Q ) J 

where r& are Gaussian distributed random numbers, <ttv is the quadratic sum of the normal- 
ization errors and at is the total error, obtained by summing in quadrature the statistical 
and systematic errors, the latter assumed to be uncorrelated. 

Following Ref. [18], the covariance matrix for experimental data points is evaluated 
using: 

COVjj = (T Nl CF Nj + Sijdij , (3.2) 

while during the fit for each replica, we minimize: 

^data / Aart),k ( n 2\ Anet),k ( n 2\\ 

x m = ^ A M) ; m 

i=i V °\,t J 

where 

of} = (l + r k n^4^ | <Ti t ■ (3.4) 

' V 4 exp) (x,Q 2 )/ ' 

The number of Monte Carlo replicas of the data is determined by requiring that the 
average over the replicas reproduces the features (central values, errors and correlations) of 
the original experimental data to a required accuracy. The quantitative check is performed 
by means of the statistical estimators described in the appendix of Ref. [2] and the results 
for sets of 10, 100 and 1000 replicas are collected in Tab. |2] for the proton target data and 
in Tab. || for the deuteron target data. We observe that all the considered estimators have 
the correct scaling behaviour as the number of replica grows. We also point out that the 
large percentage error on the deuteron central values is due to a bulk of data whose values 
are close to zero. 

3.2 Neural Networks as unbiased interpolants 

Artificial neural networks, see e.g. Ref. [19], are a class of algorithms which provide a 
robust and universal approximant to incomplete or noisy data, with the only requirement 
of continuity. Neural networks are universal approximators for measurable functions [20]. 
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10 


100 


1000 


(PE 

r 




14.20% 
0.974203 


3.21% 
0.998114 


1.83% 
0.999682 


(V [a( art )]) dat 
(PE [<> rt ))]) dat 

(^ (art) )dat 

r [a^] 


3.1 • 10" 3 

00.4070 

0.0699 
0.989956 


1.5 • 10~ 3 

10/1/10/ 

12.44% 
0.0766 
0.997808 


6.5 • 10" 4 

A onoy 
4.ZU/0 

0.0768 

0.999793 


(V [p^^dat 
(P (art) )dat 

r [p^] 


y.o • iu 
0.1469 
0.638102 


on 1 n.— 3 

o.y • iu 
0.1567 
0.944682 


no i n— 4 

y.z • IU 

0.1585 
0.993523 


(V [C0v( art )]) dat 

(cov( arf )) dat 
r [cov( ar *)] 


7.6 • 10" 5 
0.00166 
0.898803 


1.1 • 10" 5 
0.00154 
0.986219 


1.5 • 10" 6 
0.00160 
0.998858 



Table 2: Statistical estimators for Monte Carlo replicas of A\ for the proton data. The experimental 
data have (cr (oxp) } d at = 0.0764, (p( exp )) d at = 0.1566, and (cov(° xp )) dat = 0.00153. 





10 


100 


1000 


(PE 

r 


<4 art) > rcP J> 


89.71% 
0.977524 


24.90% 
0.98633 


5.97% 
0.999865 


(V [a^]) d&t 
(PE [<^ rt ))]) dat 

(^ (art) )dat 
r [<j( ar V] 


5.3 • 10" 3 
35.03% 
0.0689 

0.977501 


1.7- 10" 3 
11.75% 
0.0739 

0.997965 


6.7 • 10" 4 
4.4% 
0.0734 

0.999705 


(V [p^^dat 
(P (art) )dat 

r [p^] 


1.0 • 10" 2 

0.0878 
0.612155 


9.2 • 10" 3 

0.0904 
0.932158 


8.7 • 10" 4 

0.0861 
0.992952 


(V [cov( art )]) dat 
(cov( art )) dat 
r [cov( ar *)] 


6.3 • 10" 5 
0.00133 
0.959275 


1.0 • 10" 5 
0.00145 
0.995339 


1.0 • 10" 6 
0.00134 
0.999479 



Table 3: Statistical estimators for Monte Carlo replicas of A\ for the deuteron data. The experi- 
mental data have (cr( oxp )} da t = 0.0733, (p (oxp) ) d at = 0.0862, and (cov( oxp )) dat = 0.00135. 

This means that any continuous function can be approximated to any degree of accuracy by 
a sufficiently large neural network with one hidden layer and non-linear neuron activation 
function. 

One of the main reasons to use neural networks in place of any other redundant 
parametrization is the existence of efficient techniques for training them, i.e. determin- 
ing the parameters of the network (thresholds and weights) so that it reproduces a given 
set of input-output data. Equivalently one could say that a sufficiently large neural network 
provides a description of the data which is largely free of functional bias. 

The analysis presented here uses a class of neural networks known as multilayer feed- 



- 6 - 



forward perceptrons, trained using a genetic algorithm [21,22]. The networks we employed 
have one hidden layer and a 2-4-1 architecture, which gives us a total of 17 free parameters 
for each network to be determined during the training. The guidance principle in the choice 
of the network architecture to be used is that it should provide a redundant parametrization 
for the data to be fitted, i.e. the network should have enough flexibility to fit not only 
the underlying physical law but also the statistical fluctuations of the experimental data. 
This property is crucial in ensuring that the fit results are not biased by the specific 
parametrization. The lack of functional bias is established a posteriori by verifying that fits 
performed with networks with different architectures lead to statistically equivalent results. 
This is achieved using the statistical estimators introduced in the NNPDF Collaboration's 
studies; the results of these comparisons are presented and discussed later. 

The training of the individual 
networks to the Monte Carlo repli- 
cas is performed by minimizing the 
figure of merit given in Eq. (3.3). 
Given the extensive size and com- 
plex structure of the parameter space 
(a neural network with n param- 
eters, weights and thresholds has 
2n\ equivalent global minimum con- 
figurations), the most efficient train- 
ing algorithm turns out to be a 
genetic algorithm. The details of 
the implementation are discussed 
in Ref. [3]. 

As already pointed out in var- 
ious references the fact that we adopt a redundant parametrization and that the figure 
of merit minimized in the training procedure is monotonically decreasing might lead to 
overfitting the data: the neural network reproduces not only the underlying physical law 
but also the statistical noise of the data sample. To prevent this from happening and to 
determine the optimal fit we adopt a criterion to stop our fit based on the cross-validation 
method. Once again our procedure is completely analogous to the one used for the unpo- 
larized NNPDF fits. 

For each replica of the experimental data we subdivide the data into a training and a 
validation set, respectively containing a fraction / tr and (1 — / tr ) of randomly chosen data 
points of each experiment. 

We train one neural network on each replica of the data using the x 2 °f the training 
set as a figure of merit to be minimized. In parallel we compute the x 2 °f the validation 
set. We stop the training when we find that the x 2 smeared over a given number of 
generations is decreasing for the training set while increasing for the validation set. A 
graphical illustration of such a behaviour for one of the replicas in the reference fit is given 
in Fig. 0. 




Figure 2: \ 2 f° r the training (red) and validation (green) 
sets of one replica in the reference fit to A\. 
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4. Phenomenology 

The study of the first moments of polarized structure functions is of phenomenological 
interest, since they can be used to extract information on the fraction of polarization 
carried by partons and on physical couplings. In the MS scheme we have 



rHQ 2 



dx g± (x,Q 



(4.1) 



36 L 



a s ± 3a 3 )AC^i(a s (Q' z )) + Aa AC^(a s (Q 2 )) 



MS i 



where AC^|(a s (Q 2 )) and ACg LS (a s (Q 2 )) are the first moments of the non-singlet and 
singlet Wilson coefficient functions, respectively, and 



a 3 = (Am + Aft) - (Ad + Ad) , 

a 8 = (An + Au) + (Ad + Ad) - 2(As + As) , 

a = (An + Au) + (Ad + Ad) + (As + As) = AS . 



(4.2) 
(4.3) 
(4.4) 



Using isotopic spin invariance, it can be shown that 0,3 is the axial coupling qa = 
Ga/Gv that governs neutron /3-decay. Accurate measurements yield (see e.g. Ref. [9]): 

qa = 1.2670 ± 0.0035 . (4.5) 

The difference of the g\ moments for proton and neutron leads to the Bjorken sum 



rule 



rf (Q 2 ) = r[(Q 2 



1 



r n AQ Z ) = -gAAC^(a s (Q')) + 8 T + 5 T ; 



(4.6) 



where 6t is the target mass correction and 5t is the correction due to higher twists. Target 
mass corrections have been studied in Refs. [23-26] and can be evaluated for any moment 
n at the first order in m 2 /Q 2 using [24]: 



tT = gt\Q 2 )-9 { $(Q 2 ) 



m 2 n(n + 1) 
Q 2 (n + 2) 2 



n + 4)g^ 2 (Q 2 ) + 4'--^g^ 2 (Q 



n + 2 



n+2i 



(4.7) 



Jq dxx n 1 gi(x, Q 2 ) and is the structure function taken at zero mass 



where gi\Q 2 ) 
of the nucleon. The higher-twist contribution is simply 



(4.8) 



where U4 can be extracted from experimental data at low Q 2 such as the CLAS data [27]. 
Finally, the coefficient function of Eq. (4.6) has been calculated in Ref. [28] and up to order 
is given by: 

^2\^ 2 



AC^(a s (Q 2 )) = l-^l-(^ 



7T 



3 



a s (Q' 



TT 



115 



- 41.4399 - 7.6072 n f + n} 

1 1 648 f 



a s (Q 2 



TT 



(4.9) 
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For the running coupling we use the expanded solution of the renormalization group equa- 
tion, up to NNLO we have: 



a s (Q 2 ) = a s (Q 2 ) 



LO 



1 + a s {Q 2 ) LO [a, (Q 2 ) LO - a s (M 2 )] (b 2 - b\) 



J- rv (D 2 \ h It, US ^ 2 )nLO 



with 



(Q 2 )nlo = <** (Q 2 )lo 1 - (Q 2 ) LO In ( 1 + Poa s (Mf) In ^ 



(4.10) 



, (4.11) 



a 



(Q 2 ) 



a s (Mf) 



LO 



l + /3 « s (Mf) ln^ 
and the beta function coefficients given by 

2 



.da s (Q 2 
~dQ 2 



= -J2Pka s (Q 2 ) k+2 , a s (Q 2 ) = 



k=0 



where 



«,(Q 2 

47T 



(4.12) 



(4.13) 



38 

^ = 102 - — n/ , 

2857 5033 325 , 

ft = - Lr n / + ^ n / 



(4.14) 



and &j = f3i/f3 . 



5. Results 

In this section we present our parametrization of the proton and deuteron asymmetries 
and the structure functions extracted from them. 

We assess the quality of the fit by comparing our extraction with the experimental 
data included in the analysis, and by studying the stability of our results under variations 
of the parametrization used for the networks. 

Then, as an example of a possible application of our result to a phenomenological 
analysis, we study the extraction of the physical parameters (the strong coupling constant 
a s and the axial coupling gA) from the Bjorken sum rule. In order to give a faithful error 
on the extracted quantities we study the impact of the different assumptions which are 
needed to reconstruct the structure functions and then to evaluate the Bjorken sum rule. 
Results are compared to existing estimates. 

5.1 The final fit and its statistical features 
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Proton 


x 2 


Deuteron 


x 2 


EMC 


0.370 


COMPASS 


0.885 


SMC 


0.480 


SMC 


1.100 


SMC low-x 


1.150 


SMC low-x 


0.774 


E143 


0.904 


E143 


1.530 


E155 


0.717 


E155 


0.661 


HERMES06 


0.456 


HERMES06 


0.881 


Total 


0.666 


Total 


0.986 



Table 4: The \ 2 of the fit for proton and deuteron 
data. 



In Tab. | we show the xV-^data 
for each target and each experimental 
data set included in the present anal- 
ysis. We first observe the overall good 
quality of our fit. For the proton the 
small values of \ 2 for EMC, SMC and 
HERMES can be explained by a pos- 
sible overestimate of experimental er- 
rors. For the deuteron all the x 2 are of 
order 1, except for the E155 data set 
which has a value of \ 2 significantly 
smaller than one. The somewhat larger 

value of x 2 for the E143 deuteron data set can be understood by looking at Figs. || and [| 
where we present a comparison of our fit to experimental data in different kinematical re- 
gions. We observe that in the case of E143 the deuteron data show small incompatibilities 
among themselves, and the large value of \ 2 is a reflection of this. It is interesting to re- 
mark that a careful analysis of the x 2 value for each experiment allows the identification of 
potential incompatible data. This feature had already been pointed out in the unpolarized 
studies by the NNPDF Collaboration. 

In Tabs. ||, |(| and ^ we study the self-stability of the fit and the stability against the 
variation of the parametrization with respect to a smaller and a larger architecture. To this 
extent we define four different regions: one where we expect our fit to be an interpolation 
of the available data (Data region) and three where its behaviour is extrapolated to regions 
of the (x, Q 2 )-plane not covered by present data: 

• Data: 0.01 < x < 0.75 and 2 GeV 2 < Q 2 < 20 GeV 2 ; 

• Low-x: 0.0001 < x < 0.001 and 2 GeV 2 < Q 2 < 20 GeV 2 ; 

• Low-Q 2 : 0.2 < x < 0.8 and 0.1 GeV 2 < Q 2 < 2 GeV 2 ; 

• High-Q 2 : 0.2 < x < 0.8 and 20 GeV 2 < Q 2 < 60 GeV 2 . 

We observe that all the estimators for self-stabilities are of order unity (or smaller), mean- 
ing that different subsets within the whole ensemble of replicas have the same statistical 
features. 

When we compare our final fit to a fit performed using networks with a smaller ar- 
chitecture, we notice that the the two fits are statistically equivalent. The same happens 
for the comparison with a fit done with networks with a larger architecture, with the only 
exception of the errors on the deuteron fit in the extrapolation (all distances are order 1.5), 
which show some minor instability. 

5.2 Structure functions reconstruction 

In order to reconstruct the structure function g\ from data on the asymmetry A\ as given 
in Eq. fl2,4|) some additional assumptions are needed. In the following we assess the impact 
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Figure 3: The fitted asymmetries compared to proton (left) and deuteron (right) data for 



0.01 GeV 2 < 



< 1 GeV 2 (upper row), 1 GeV 2 < 



< 3 GeV (central row) and 



3 GeV 2 < Q 2 < 5 GeV 2 (lower row). In the plots A\ is evaluated at the central value of each 
Q 2 range. 



of our assumptions for g2 , F% and R on the determination of first moment of g\. These 
checks are done using an ensemble of 100 replicas, which is enough to this purpose, and in 
a range of x and Q 2 which is entirely in the data region in order to avoid any extrapolation 
effects. Finally, we compare our result for 1000 replicas with the sum rules obtained by 
experimental collaborations. 

The first assumption whose impact we consider is the one on the structure function 
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Figure 4: The fit compared to proton (left) and deuteron (right) data for 5 GeV 2 < Q 2 < 10 GeV 2 
(upper row), 10 GeV 2 < Q 2 < 30 GeV 2 (central row) and 30 GeV 2 < Q 2 < 60 GeV 2 (lower row). 
In the plots A\ is evaluated at the central value of each Q 2 range. 



#2, which is evaluated from the Wandzura-Wilczek relation [29] 

9 Y w (x,Q 2 ) = - gi (x,Q")+ ( X % x {y,Q 2 ). 

Jx y 

Inserting this expression into Eq. (|2.4| ) gives 

9i(x, Q 2 ) = (a^x, Q^F^x, Q 2 ) + 7 2 jf * %i(y, Q 2 )) , 



(5.1) 



(5.2) 
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Proton 


Data 


Low-x 


Low-Q 2 


High-Q 2 


(d[Ai]) 
(d[<r Al ]) 


0.963 ±0.011 
0.840 ± 0.006 


0.660 ±0.011 
0.618 ±0.011 


1.085 ±0.015 
0.966 ±0.014 


0.966 ±0.014 
0.905 ±0.012 


Deuteron 


Data 


Low-x 


Low-Q 2 


High-Q 2 


(d[A,}) 
(d[a Al ]) 


0.772 ± 0.008 
0.818 ±0.007 


0.670 ± 0.012 
0.899 ±0.014 


0.804 ±0.012 
0.730 ±0.011 


0.688 ±0.011 
0.773 ±0.011 



Table 5: Self-stability estimators evaluated with 100 replicas. The entries in the table show the 
statistical differences between results based on different subsets of 100 replicas randomly chosen in 
our Monte Carlo ensemble. 



Proton 


Data 


Low-x 


Low-Q 2 


High-Q 2 


(d[Ai]) 
(d[*Ai]) 


0.952 ±0.010 
1.104 ±0.008 


0.792 ±0.014 
1.405 ±0.013 


0.859 ±0.012 
0.975 ±0.014 


1.295 ±0.016 
1.002 ±0.011 


Deuteron 


Data 


Low-x 


Low-Q 2 


High-Q 2 


(d[a Al }) 


1.217 ±0.012 
0.963 ± 0.008 


1.302 ±0.012 
1.199 ±0.010 


1.324 ±0.017 
0.962 ±0.010 


0.703 ± 0.010 
1.689 ±0.019 



Table 6: Stability estimators for the reference fit (architecture 2-4-1) compared to a fit with a 
smaller architecture (2-3-1). 



Proton 


Data 


Low-x 


Low-Q 2 


High-Q 2 


(d[M]) 
(d[<r Al ]) 


1.076 ±0.012 
1.258 ±0.010 


1.179 ±0.017 
1.354 ±0.014 


0.693 ±0.011 
0.709 ±0.010 


1.625 ±0.018 
1.151 ±0.016 


Deuteron 


Data 


Low-x 


Low-Q 2 


High-Q 2 


(d[Ai]) 
(dM) 


0.867 ± 0.009 
1.460 ± 0.009 


0.794 ± 0.012 
1.941 ±0.015 


0.856 ±0.015 
1.410 ±0.011 


0.840 ±0.017 
1.303 ±0.011 



Table 7: Stability estimators for the reference fit (architecture 2-4-1) compared to a fit with a 
larger architecture (2-5-1). 



which needs to be evaluated iteratively. To this purpose we take the initial value gi(x, Q 2 ) 
evaluated with g2(x, Q 2 ) = 0, and we use 

$™\x, Q 2 ) = ^ U(x, Q 2 )Fi(x, Q 2 ) + 7 2 j l %?™-V {y , Q 2 )\ . (5.3) 

From Tab. || we see that for Q 2 values above 2 GeV 2 one iteration is enough to stabilize 
the result for the first moment of g\ computed in the data region. For lower scales, say 
Q 2 ~ 1 GeV 2 , at least two iterations of the Wandzura-Wilczek relation are needed in order 
to obtain a stable result. In the following the index iww will be omitted as the number 
of iterations used should be evident from the scale at which the first moment of g\ is 
evaluated. 

For the unpolarized structure function F% we use the parametrization given in Ref. [2] 
for the proton and the one given in Ref. [1] for the deuteron. Since these parametrizations 
have also been extracted using a Monte Carlo procedure, ensembles of replicas are available 
for F2; hence the result for g± is evaluated as: 



JV rcp 

rep fc=i 



(5.4) 
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Proton 


Q 2 = 2 GeV 2 


Q 2 = 5 GeV 2 


Q 2 = 10 GeV 2 


Q 2 = 20 GeV 2 


iww — 
iww — 1 

iww 2 


0.1184 ± 0.0069 
0.1086 ±0.0062 
0.1072 ± 0.0061 


0.1196 ± 0.0068 
0.1154 ±0.0066 
0.1152 ± 0.0065 


0.1272 ± 0.0106 
0.1251 ±0.0104 
0.1251 ±0.0104 


0.1401 ± 0.0216 
0.1391 ±0.0215 
0.1391 ±0.0215 


Deuteron 


Q 2 = 2 GeV 2 


Q 2 = 5 GeV 2 


Q 2 = 10 GeV 2 


Q 2 = 20 GeV 2 


iww = 
iww = 1 
iww = 2 


0.0459 ± 0.0049 
0.0414 ± 0.0044 
0.0407 ± 0.0042 


0.0414 ± 0.0036 
0.0396 ± 0.0035 
0.0395 ± 0.0034 


0.0401 ± 0.0064 
0.0392 ± 0.0064 
0.0392 ± 0.0064 


0.0397 ± 0.0149 
0.0393 ± 0.0149 
0.0393 ± 0.0149 



Table 8: First moment for x between 0.01 and 0.75 for different number of iterations of the 
Wandzura-Wilczek relation. 



which takes into account both the uncertainty on A\ and the one on F2 (with g^ {x,Q 2 ) 
we denote the expression in Eq. ( |5,1| ) evaluated for the k-th replica). Since there is no 
correlation between the extraction of A\ and the one of F% the replicas of A\, and F2 can 
be sampled independently. 

In order to estimate the contribution of the uncertainty on F2 to the uncertainty on 
<7i, we can recompute g\ as 



gi(x,Q 2 



1 



rep 



N Iep 

£ 

k=l 



A k) (x,Q 2 



(1+7 2 ) 



2x[l + R(x,Q 2 )] 



(F2)(x,Q 2 ) + 1 2 gi k \x,Q 



(5.5) 



where for each k-th. replica of A\ we use the averaged value of the unpolarized structure 
function 



iVrcp 

{F2)(x,Q 2 ) = -^Y. F 2 k) ( x >Q 2 y 

JYrop k=l 



(5.6) 



This procedure clearly freezes the fluctuations in F2, which is kept fixed to its average 
value. The result is given in Tab. ||, where we see that the contribution to the uncertainty 
on the first moment of g\ due to F2 is negligible. In the following we will always use g\ as 
given from Eq. (|5.4|). 



Proton 


Q 2 = 2 GeV 2 


Q 2 = 5 GeV 2 


Q 2 = 20 GeV 2 


Eq. ( 
Eq. ( 


5.4 
5.5 


) 
) 


0.1086 ± 0.0062 
0.1086 ±0.0059 


0.1154 ±0.0066 
0.1154 ±0.0064 


0.1391 ±0.0215 
0.1391 ±0.0213 


Deuteron 


Q 2 = 2 GeV 2 


Q 2 = 5 GeV 2 


Q 2 = 20 GeV 2 


Eq. ( 
Eq. ( 


5.4 
5.5 


) 
) 


0.0414 ± 0.0044 
0.0414 ± 0.0038 


0.0396 ± 0.0035 
0.0396 ± 0.0034 


0.0393 ± 0.0149 
0.0394 ±0.0150 



Table 9: First moment for x between 0.01 and 0.75 with and without the error on F2 

Finally a parametrization of R(x,Q 2 ) is needed in order to extract g\ from A\. Here 
we use Rslac(x,Q 2 ) given in Ref. [30,31]. Such a parametrization provides also an error 
estimate, which we use to assess the impact of Rslac(x, Q 2 ) on the total uncertainty of the 



first moment of g\. In Tab. 10 we compare the sum rule evaluated with the central value 
of Rslac(x, Q 2 ) with the one obtained by taking into account the error on Rslac(x, Q 2 )- 
This is achieved by letting Rslac(x, Q 2 ) fluctuate within its own error in the Monte Carlo 
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sample; for the k-th replica we use: 



Rslac(x, Q 2 ) + r^AR S LAc(x, Q 2 ) , 



(5.7) 



where ARslac( x > Q 2 ) is the error on the parametrization, and A k ^ is a univariate Gaussian 
random number. Since Rslac{x,Q 2 ) is a parametrization of experimental data, we take 
the error as a statistical one, with no correlation between different replicas, and thus we 
use a different random number each time a value of Rslac(x,Q 2 ) is needed. From the 
results collected in Tab. 10 we conclude that the error on Rslac(x,Q 2 ) is also negligible. 



Proton 


Q 2 = 2 GeV 2 


Q 2 = 5 GeV 2 


Q 2 = 20 GeV 2 


Rslac(x, Q 2 ) 

Rslac(x,Q 2 ) + r {k) ARslac(x,Q 2 ) 


0.1086 ± 0.0062 
0.1086 ± 0.0062 


0.1154 ±0.0066 
0.1154 ± 0.0066 


0.1391 ±0.0215 
0.1392 ± 0.0216 


Deuteron 


Q 2 = 2 GeV 2 


Q 2 = 5 GeV 2 


Q 2 = 20 GeV 2 


Rslac(x,Q 2 ) 

Rslac(x,Q 2 ) + r (fc) ARslac(x,Q 2 ) 


0.0414 ± 0.0044 
0.0415 ± 0.0044 


0.0396 ± 0.0035 
0.0397 ± 0.0035 


0.0393 ± 0.0149 
0.0395 ± 0.0148 



Table 10: First moment for x between 0.01 and 0.75 with and without the error on R 



We will now compare our results for the integral of g\ at different scales and over 
different x ranges obtained using our ensemble of 1000 replicas with those obtained by 
different experimental collaborations. 



In Tab. 11 results for the proton and 
the deuteron sum rules are compared to the 
result of Ref. [11]. We observe that the 
results are compatible within errors, and 
that our evaluation has a larger error. 



Target 


SMC98 


This Analysis 


Q 2 = 10GeV 2 


P 


0.131 ±0.009 


0.139 ±0.015 


d 


0.037 ±0.007 


0.035 ±0.011 



Table 11: Comparison of the proton and 
deuteron sum rules ( Jqq 03 dxgi(x, Q 2 )j as deter- 
mined in the present analysis with the results ob- 
tained by the SMC collaboration [11]. 



In Tab. 12 we compare our result with 
the ones in Refs. [13], and [15]. First we 
notice that the errors are of the same size, 
while our central values are systematically 
smaller, with a significant difference for the 

proton at low Q 2 . A substantial part of effect can be attributed to the different parametriza- 
tion used for the unpolarized structure function. Indeed, if we evaluate the sum rule of E143 
with the SMC98 F% parametrization [11,32,33], at Q 2 = 2 GeV 2 we obtain 0.116 ± 0.008 
which is less than one sigma away from the result in Refs. [13]; the same happens for 
the HERMES06 case, using the ALLM parametrization [34] at Q 2 = 2.5 GeV 2 we get 
0.1188 ± 0.0073. This can be understood looking at Fig. |5| where we compare the different 
parametrizations for F% used in the different analysis. 

It is clear that, while the different F% parametrizations agree in the kinematical region 
covered by experimental data, they differ significantly at low-Q 2 in the large- x region 
where there are no data and an extrapolation is needed. For the ALLM and the SMC98 
parametrizations the large- x behaviour is determined by the chosen functional form; the 
NNPDF parametrization interpolates by continuity from the last experimental point to 
the kinematical constrain F2(x = 1, Q 2 ) = 0. The difference among the parametrizations 
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f;<x,q 2 ) |f;<x,q 2 ) 




X 



Figure 5: Comparison of NNPDF, SMC98 and ALLM parametrizations of the unpolarized struc- 
ture function F% in the region where we evaluate the Bjorken sum rule. 



is then enhanced once we multiply by the asymmetry A\ to reconstruct the polarized 
structure function g\. 



Target 


E143 


This Analysis 


Q 2 = 2GeV 2 


P 


0.120 ±0.007 


0.102 ±0.007 


d 


0.047 ±0.006 


0.042 ± 0.005 


n 


-0.022 ±0.013 


-0.011 ±0.011 


NS 


0.149 ±0.016 


0.113 ±0.016 


Q 2 = 5GeV 2 


P 


0.116 ±0.007 


0.106 ± 0.006 


d 


0.043 ±0.004 


0.040 ± 0.003 


n 


-0.025 ±0.009 


-0.018 ±0.009 


NS 


0.141 ±0.013 


0.124 ±0.014 



Target 


HERMES 


This Analysis 


Q 2 = 2.5GeV 2 


P 


0.1201 ±0.0090 


0.1055 ± 0.0066 


d 


0.0428 ± 0.0035 


0.0416 ± 0.0043 


n 


-0.0276 ± 0.0093 


-0.0154 ±0.0107 


NS 


0.1477 ±0.0167 


0.1209 ±0.0152 


Q 2 = 5GeV 2 


P 


0.1211 ±0.0092 


0.1097 ±0.0065 


d 


0.0436 ± 0.0035 


0.0407 ± 0.0033 


n 


-0.0268 ± 0.0094 


-0.0218 ±0.0093 


NS 


0.1479 ±0.0169 


0.1315 ±0.0144 



Table 12: Comparison of the integral of g± over different x ranges, at different scales, as determined 
form the present analysis with the results obtained form by E143, left pad: J 03 dxg\{x, Q 2 ) 
HERMES, right pad: J ° 9 21 dxgi(x, Q 2 ) 



and 



Since no neutron target data have been used in our fit, the neutron structure function, 
<7™, is evaluated from the proton and deuteron ones as 

gUx,Q 2 ) = 2 f} X ^ ) D -<f 1 (x,Q 2 ), (5.8) 

where ujd is the probability of the deuteron to be in the D state; we use u>d = 0.05 which 
covers most of the published values [35]. We observe that, even if the neutron sum rule is 
a pure prediction, it is compatible with other estimations. 

In Fig. ^ we show a comparison of the polarized structure function as extracted in this 
analysis to data, in the region where we evaluate the Bjorken sum rule. The comparison 
shows a good agreement. 
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Figure 6: Plot of the structure functions in the region where we evaluate the Bjorken sum rule. 
The fit curves are taken at Q 2 = 5 GeV 2 . 

5.3 Extraction of couplings 

In order to extract the strong coupling a s and the axial coupling g^ from the values of the 
Bjorken sum rule we need to extrapolate our fit in the Bjorken variable x down to x = 
and up to x = 1. 

In this section we discuss the impact of these extrapolations on the extraction of the 
couplings, and we assess the impact of target mass corrections. Finally we present the 
results we obtain for a s and gA from our fit. All checks are performed with 100 replicas, 
while for final results we use the full set of 1000 replicas. 

The extrapolation at large-x is embedded in the parametrization of as discussed 
in the previous section. Therefore we do not need any further assumption to constrain the 
large-x behaviour. 

The low— x behavior of the structure function g% is instead very weakly constrained by 
data, and the Regge behaviour is usually assumed; following Ref. [36] we write: 



gi(x,Q 2 ) ~Ax b 



(5.9) 



with < b < 0.5. Such an assumption requires to choose a value of x matc h such that for 



x < a^match the Regge behaviour is assumed to set in. The normalization factor A in Eq. 5.9 
is then determined by the matching condition 



9i 



(fit), 



^matchj Q ) — A X matc j 1 . 



(5.10) 



In order to choose the matching point, we fix the Regge exponent to 0.2 and then we 
proceed as in Ref. [37]: we evaluate the integrals at different values of Q 2 in the range 
< x < 1 for different choices of £ ma tch> and we look for the minimum value of the error 
for each value of Q 2 . 



From the results collected in Tab. 13 we see that x ma t c h grows as Q gets larger. This 



is understood looking at Fig. |l] where we see that for larger scales the coverage of the data 
moves towards higher values of x. 
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Once the matching point is been de- 
termined, in order to take into account 
the uncertainty on the value of the Regge 
exponent and the one on the choice of 
the matching point, we randomize the 
Regge exponent in the range —0.1 < b < 
0.5 and we choose the matching point to 
be in the range x matc h < x < 2x ma , tc h- 



In Tab. |14j we present the compari- 
son for the first moment of the structure 
function g\ evaluated with and without 
the target mass correction as given in 
Eq. (|4.7|) . We observe that the shift on 
the values of the moment due to the in- 
clusion of these effects is smaller than 
the experimental error even at the low- 
est Q 2 . 

In principle we could extract g^, a s 



Q 2 (GeV 2 ) 


■^match 


rNs (Q2) 


Error 


1 


0.0100 


0.12499 


0.020989 


2 


0.0100 


0.1356 


0.018239 


3 


0.0100 


0.14324 


0.017827 


4 


0.0200 


0.13847 


0.018275 


5 


0.0200 


0.14322 


0.019021 


6 


0.0200 


0.14757 


0.020115 


7 


0.0200 


0.15142 


0.021429 


8 


0.0300 


0.1458 


0.02262 


9 


0.0300 


0.14827 


0.023859 


10 


0.0300 


0.15054 


0.025165 


11 


0.0300 


0.15265 


0.02653 


12 


0.0500 


0.14111 


0.027703 


13 


0.0500 


0.14241 


0.028749 


14 


0.0500 


0.14366 


0.029828 


15 


0.0500 


0.14487 


0.030941 


16 


0.0500 


0.14604 


0.032087 


17 


0.0500 


0.14718 


0.033267 


18 


0.0500 


0.14829 


0.03448 


19 


0.0800 


0.13397 


0.03552 


20 


0.0800 


0.13461 


0.036394 



Table 13: First moment of g\ with different values 
of Q 2 : the error on T^ S (Q 2 ) has been added a 100% 
uncertainty on the low- a; extrapolation. 



and the higher-twist term by fitting Eq. (|4.6 
evaluated from data at a given value of 
Q 2 . In practice we evaluate Nq2 dif- 
ferent moments taken at different Q 2 in 
the kinematical region where we have 
a good coverage by experimental data: 
2 GeV 2 < Q 2 < 20 GeV 2 . Indeed for 
Q 2 > 20 GeV 2 the errors on the com- 
puted moments become so large that their weighted contribution in the combination is 
negligible. On the lower side of the energy range we choose to start from Q 2 = 2 GeV 2 , 
since below this scale a perturbative QCD approach might not be reliable. For this reason 
we do not fit the higher-twist term, but we will access its contribution by varying the lower 
cut in Q 2 . 

We then proceed following the 
procedure described in detail in Sect. 
4.3 of Ref. [38]: the extraction of 
couplings is done by combining mo- 
ments at different values of Q 2 in the 
chosen range and fitting Eq. ( |4,6j ) 
using MINUIT [39] where g A and 
a s are the chosen as free parame- 
ters. The moments at different Q 2 
are correlated, since they are com- 



Proton 


Q 2 = 1 GeV 2 


Q 2 = 2 GeV 2 


noTMC 
wTMC 


0.1410 ±0.0109 
0.1459 ±0.0116 


0.1264 ± 0.0083 
0.1276 ± 0.0086 


Deuteron 


Q 2 = 1 GeV 2 


Q 2 = 2 GeV 2 


noTMC 
wTMC 


0.0493 ±0.0141 
0.0515 ±0.0158 


0.0350 ± 0.0058 
0.0355 ± 0.0059 



Table 14: First moment with and without TMC with 
b = 0.2 for the low- a; extrapolation matched at x = 0.001. 



puted using the same fitted parametrization. As detailed in Ref. [38] these correlations 
induce numerical instabilities in the inversion of the correlation matrix and off-diagonal 
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instabilities due to non-diagonal elements in the correlation matrix becoming dominant. 
Both these instabilities lead to unreliable results for the extracted couplings. 

In order to fix the maximum value of Nq2 
for which the extraction of the parameters is 
numerically stable and reliable, we study the 
error on the determination of g A and a s as we 
vary the number of included moments. Once 
we exclude moments with large correlations, we 
are left with a small number of combinations, 
which are showed in Tab. 

The combination giving the smallest error 
(Q 2 = 2, 5 GeV 2 ), once we evaluate asymmetric 
errors, yields 



Q 2 


9A 


a s {M%) 


2+5 


1.04 ± 0.12 


0.126 ± 0.005 


2+6 


1.07 ± 0.16 


0.127 ± 0.007 


2+7 


1.08 ± 0.20 


0.128 ± 0.008 


2+8 


1.02 ± 0.21 


0.125 ± 0.012 


3+6 


1.09 ± 0.13 


0.131 ± 0.005 


3+7 


1.12 ± 0.17 


0.132 ± 0.007 


3+8 


1.04 ± 0.20 


0.127 ± 0.013 


4+8 


1.03 ± 0.17 


0.127 ± 0.015 



0.126 



+0.004 
-0.009 ' 



(5.11) 



Table 15: Fits with different choices of Q 2 



for the strong coupling, while the error on g A is found to be symmetric. 

The only sources of theoretical uncertainty 
left to consider are the one due to the choice of 
factorization scale Q = kpm q , which we study 
by varying kp in the range 0.5 < kp < 2 and 
to the higher-twist contribution. The results 
of the variation of the factorization scale are 
To take into account the 



9A 


a s (M 2 ) 


kp 


1.02 ± 0.12 


0.126 ± 0.006 


0.5 


1.04 ± 0.12 


0.126 ± 0.005 


0.75 


1.01 ± 0.11 


0.121 ± 0.006 


1.5 


1.03 ± 0.12 


0.120 ± 0.005 


2.0 



shown in Tab. 16. 



Table 16: Reference fit (Q 2 = 2,5 GeV 2 , 
NNLO, 1000 reps, kp — 1) compared with 
variations of the factorization scale. 



higher-twist contribution we take as an esti- 
mate the variation of the central values once 
the lower Q 2 value is moved up to Q 2 = 3 GeV 2 
(see Tab .|15|) and down to Q 2 = 1 GeV 2 (g A = 
0.99 ± 0.14 and a s (M 2 ) = 0.117 ± 0.007). 

In conclusion, we obtain the following result for the determination of the axial coupling 
gA and the strong coupling constant a s 



g A = 1.04±0.12(exp.)+0;^(theo.) 
a s (M 2 ) = 0.126l8;^(exp.)18:8?f (theo.) 



1.04±0.13(tot.) 

= o.i26iHSl(tot.; 



(5.12) 



which are compatible with previous extractions [42,43] from polarized DIS and the Bjorken 
sum rule. 



6. Conclusions 

We extracted a parametrization of the spin asymmetries A\ ,d , based on all available DIS 
data using the Monte-Carlo sampling techniques and neural networks as basic interpolation 
tools. We checked in the process that the statistical methods developed for the unpolarized 
studies by the NNPDF Collaboration can be naturally extended to handle the new data 
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sets considered in this work. Our main result is an effective tool, which we used to test 
different assumptions needed to reconstruct the polarized structure function g± . As an 
example of possible applications we compared to previous estimations of experimental sum 
rules, and we found that the used parametrization for the unpolarized structure function 
F2 can be a sizable source of error at low values of Q 2 . We also performed a study of 
the Bjorken sum rule, and the extraction of the axial coupling and the strong coupling, 
obtaining values which are compatible with previous analysis. 

It would be interesting to compare the results obtained for the Bjorken sum rule when 
determined from global QCD fits to polarized DIS, SDIS and hadron-hadron collisions 
data, like the one presented in [44], especially once W production data from RHIC will be 
included in such fits, providing an extra constraints on the light flavours separation. 

The present study is also meant to be a first step towards the application of the NNPDF 
techniques to the determination of a set of polarized PDFs with a faithful error estimation. 
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A. Experimental errors 

Experimentally, we have 

C N--N. 

A\\ ~ 

where 

• C is a nuclear correction that depends on the material the target is made of; 

• / is the dilution factor which accounts for the fact that only a fraction of the target 
nucleons is polarizable; 

• Pf, and Pt are the beam and target polarizations; 

• N_r + \ is the number of scattered electrons/muons per incident charge for negative 
(positive) beam helicity. 

Thus, most of the errors quoted by experiments are normalization errors. 

• EMC [10]: A2 = is assumed; 9.6% overall normalization due to beam and tar- 
get polarization; multiplicative errors on R and /; additive errors on A2, the false 
asymmetry K and the radiative correction. 
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(A.l) 



• SMC98 [11]: A2 = is assumed; multiplicative errors on Pt, Pb, R, f and the 
polarized background APb g ; additive errors on A2, the false asymmetry AAf a i se , the 
radiative correction and the momentum resolution. 

• SMC low-x [12]: A2 = is assumed; multiplicative errors on Pt, Pb, R, f and the 
polarized background APb g ; additive errors on A2, the false asymmetry AAf a i se and 
the radiative correction. 

• E143 [13]: 52 is evaluated using the Wandzura-Wilczeck relation 



using and empirical fit of g±/Fi = ax a (l + bx + cx 2 )(l + Cf(Q 2 )); multiplicative 
errors on P t , Pb, f and the nuclear correction C which account for a total 3.7% for 
the proton and 4.9% for the deuteron; additive uncorrelated error on the radiative 
corrections. 

• E155 [14]: 02 is evaluated in the same way of E143, but the parameters of the fitted 
functional form have different values; we will add as a shift the difference between 
A\ and gi/F±; multiplicative errors on Pt, Pf,, f and the nuclear correction C which 
account for a total 7.6% for the proton; additive uncorrelated error on the radiative 
corrections. 

• HERMES06 [15]: a parametrization for §2 is fitted to existing data; normalizations 
errors of 5.2% for the proton and 5% included in the systematic quoted for each data 
point; additional additive error on the parametrization used for $2- 

• COMPASS [16]: A2 = is assumed; multiplicative errors on P t , Pb, the dilution 
factor / and the depolarization factor D; additive errors on the false asymmetry and 
the radiative correction. 
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